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Abstract 

We extend a path-integral approach to bosonization previously developed in 
the framework of equilibrium Quantum Field Theories, to the case in which time- 
dependent interactions are taken into account. In particular we consider a non 
covariant version of the Thirring model in the presence of a dynamic barrier at 
zero temperature. By using the Closed Time Path (Schwinger-Keldysh) formalism, 
we compute the Green's function and the Total Energy Density of the system. 
Since our model contains the Tomonaga Luttinger model as a particular case, we 
make contact with recent results on non-equilibrium electronic systems. 
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1 Introduction 



The study of a quantum field theory involving time- dependent perturbations is a relevant 
issue for a large variety of phenomena such as formation and growth of domains in high 
energy physics pQ , decoherence and photon pair creation in a homogeneous electric field 
0, electron transport in solids 3\, etc. In recent years there has been an intense activity 
focused on the analysis of non equilibrium situations in the context of low dimensional 
electronic systems jl]. In particular, theoretical work on one dimensional (ID) structures 
relies heavily on field theories obtained as continuum limits of original lattices and on the 
systematic use of bosonization and renormalization group methods [S]. In this context, 
an alternative path-integral bosonization method, valid for the description of systems 
in equilibrium, was first established in [SI- In view of the interesting open problems in 
ID non-equilibrium fermionic systems, in this work we generalize the above mentioned 
functional bosonization method to account for the time evolution of the system. To this 
end we consider a generalized (non covariant) Thirring like model in the presence of a dy- 
namic barrier at zero temperature. By using the Closed Time Path (Schwinger-Keldysh) 
formalism (CTP)[7J, we compute the Green's function and the Total Energy Density 
(TED) of the system. Our model contains the Tomonaga-Luttinger model (paradigm of 
Luttinger liquid behavior) as a particular case [Hj. This allows us to make contact with 
recent interesting results on non-equilibrium many-fermion electronic systems fT] which 
were obtained through operational bosonization. These authors were able to obtained 
the TED for a generic impurity (without specifying its spatial dependence) by fixing 
from the beginning a spatial point of interest and disregarding the adiabatic switching 
of the time-dependent interaction. Although this procedure is fully justified on physical 
grounds, taking into account the temporal step function brings additional difficulties to 
the calculations and it is not evident a priori that results should coincide, unless a large 
time limit is considered. Then, as a secondary task, we have taken into account the 
influence of the adiabatic switching, at least at the level of the Green's function. This, 
in turn, prevented us from obtaining a TED for an arbitrary impurity geometry. We 
were then forced to specify a particular barrier, but, on the other hand, we could get the 
spatial dependence of the energy distribution. 

The paper is organized as follows. In Section 2 we show how to conciliate the path- 
integral approach to bosonization based on decoupling changes in the functional measure 
with the procedures of the Closed Time Path formulation. In Section 3 we compute the 
Green's functions and specialize our results to the case of an adiabatically switched bar- 
rier of width a and height A which oscillates with frequency Q. In Section 4 we compute 
the Total Energy Distribution function and its temporal average, in the large distance 
(time) approximation. This enables us to analyze the behavior of this distribution as 
function of energy, strength of the dynamical barrier and of the electron-electron inter- 
actions. We also comment on the spatial distribution for given energies. Our results are 
in qualitative agreement with the findings of ref. [O], which were obtained by using con- 
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ventional (operational) bosonization. In Section 5 we summarize our results and gather 
our conclusions. 

2 CTP Formalism and Functional Bosonization 

In the S matrix formalism the time is taken over the interval (— oo, +00). The state at 
t = —00 is well defined as the ground state of the non-interacting system. The interactions 
are turned on slowly. At t ~ the fully interacting ground state is -?/>(0) = S(0, —00) <3>o- 
In condensed matter physics the state at t — > 00 must be defined carefully. If the 
interactions remain on, then this state is not well described by the non-interacting ground 
state $o- Alternately, one could require that the interaction turn off at large times, which 
returns the system to the ground state <3>o- As it is well known, in 1961 Schwinger [7j 
suggested another method of handling the asymptotic limit t — > 00. He proposed that 
the time integral in the S matrix could be separated in two pieces: one going from —00 
to a certain time r) and the second one going from r to — 00 (eventually one could 
also let r — > 00). Thus the integration path is a time loop C, which starts and ends 
at t = —00, and it is composed by two branches C + (— 00, r) and C_ (r, — 00). The 
advantage of this method is that one starts and ends the S matrix expansion with a 
known state ip{— 00) = $0, which usually is the only ground state one knows exactly. For 
equilibrium phenomena the time loop method of evaluating the S matrix gives results 
that are identical to the other methods. Perhaps its main practical advantage resides in 
the possibility of describing non-equilibrium phenomena using Green's functions. The 
equation of motion for the Green's function can be cast into the form of a quantum 
Boltzmann equation for the transport theory [Hj. A disadvantage of the time- loop method 
is that it employs four different Green's functions. They are all correlation functions 
which relate the field operator of the particle at one point x = (x, t) in space-time to the 
conjugate field operator at another point x' = (x',t'): 

i(&(x!)*(yL)), teC_, t'eC + 

-i(*(x)*V)) 5 t e C+, t' e C_ 
-i(T*(x)& (*!)), t,t' eC- 

-z^f^x^V)) , M' eC + (2.1) 

where T and T denote the time and anti-time ordering operations respectively, and C + 
and C_ are the time integration paths mentioned above, i.e. they are the two branches of 
the loop C. The subscripts for the Green's functions denote the branches of the contour 
in the complex time plane on which the appropriate time coordinates reside. 

After this sketch of the CTP method, our first task is to show how this powerful 
formulation can be incorporated in the path-integral approach to bosonization based on 



GL+(x,x') = 

G+_(x,x') = 

GL_(x,x') = 

G'++(x,x / ) = 
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decoupling changes in the functional integration measure |U|. To this end we consider a 
non covariant version of the Thirring model which has been previously used to describe 
the low-energy physics of a system of non relativistic ID electrons. In order to illustrate 
the way in which a non equilibrium situation can be managed in our approach, we add 
an interaction that depends explicitly on time. This problem gives us the opportunity 
to go beyond a purely academic computation by choosing a type of perturbation that 
involves interesting physical effects. This is the case of dynamic barriers that can be also 
viewed as time-dependent spatially localized impurities [TT] [T2] . 
We start from the Euclidean action given by 

S = Sq + Sm + Si, (2.2) 

where So is the unperturbed action (in the condensed matter context it is thought of as a 
linearized free dispersion relation valid to explore the low energy-long distance physics): 

S = I d 2 x^i$^, (2.3) 
Jc 

S int describes the forward scattering of spinless fermions (electrons): 

Sint = -\\/* A(^)(x) V w (x,y) (% M *)(y) (2.4) 

in which the fermionic currents ^7^^/ are coupled through bilocal distance-dependent 
potentials, VL\(x, y). In terms of these potentials one can make direct contact with 
the forward-scattering sector of the "g-ology" model currently used to describe different 
scattering processes characterized by coupling functions g%, gi, g% and g± Neglecting 
processes associated to large momentum transfers, only the forward scattering couplings 
gi and g± play a role. The relation between these strengths and our potentials are given 
(in Fourier space) by 

<fe(p) =|(V(i)(p) + ^ ( o)(p)) 

9a(p) =f(V (0 )(p)-^(i)(p)). (2.5) 

Finally, Si describes the interaction between the electrons and a localized time- 
dependent perturbation: 

Si= I rf 2 x^ 7o A(x,t)*, (2.6) 
Jc 

where 70 is a Dirac matrix and A(x,t) contains the details of the perturbation, i.e. the 
way in which the interaction is switched on in time and the form in which it is localized in 
space. We will make a particular choice for this function later. In all cases, J c indicates 
that the integration is defined along the time contour that we briefly described above. 
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Let us now briefly show that all the steps followed in j^J, in order to express the 
functional integral associated to (|2.2|) as a fermionic determinant, are also valid when the 
time-dependent perturbation is taken into account. We consider the vacuum to vacuum 
functional: 

Z = I V^Vme~ s (2.7) 
Jc 

and recall that the electron-electron interaction piece of the action can be written as 

1 r * 



S int = -- I d'xJ^K^ (2.8) 



c 



where J M is the usual fermionic current, and is a new current defined as 

K,= f A^)(x,y)J M (y) (2.9) 
Jc 

Using a functional delta and introducing auxiliary vector fields and in the path 
integral representation of Z (See jH] for details) one obtains: 

Z = N J VVWrDApVBp exp j- J d 2 x [^{0 + 7o A (x, t))V + A^B^ + -^(A^ + J 

(2.10) 



If we define 



with satisfying 



fl M (x) = f d 2 yV,(y, X )B,(x) (2.11) 
Jc 

B M (x) = / rf , «V- 1 (y l x)S M (x) (2.12) 
Jc 

d 2 yV-\-K iy )V^y) = S(x - z)S c (t x -t z ), (2.13) 



c 



where 5c is the delta function extended to the C contour as: 



5(t-t') t,t' on C+ 
5 c {t-t')={ -6{t-t') t,t' on C_ (2.14) 
otherwise 



and changing auxiliary variables in the form 

1 
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A^ = —(A^ + BJ (2.15) 



B, = -^(i„-J3 M ) (2.16) 



we get 



Z=M J V^V^VA^VB^exp-^'^exp^- J d 2 x + 7 A(x, t) - y4)*]J (2.17) 



where 



= \\ d 2 xd 2 yV-\x,y){A,( X )AM - B,( X )B,(y)] (2.18) 

The jacobian associated with the change (A, B) — > (A, B) is field-independent and 
can then be absorbed in the renormalization constant M . Note also that the field B^ is 
completely decoupled from both the field A^ and the fermion field and then its contribu- 
tion (which describes a negative metric state) can be factorized and absorbed in M . We 
thus find 

Z=M J VA^exp(-S[A^\) det(^-4 + 7o A) (2.19) 



where 

S[A,} = \J c d 2 x d 2 y A^)V { -]{^y)A,{y) (2.20) 

Having expressed Z in terms of a fermionic determinant we are ready to apply the 
machinery of the decoupling approach to functional bosonization which is based on the 
following transformation 

tf( x ) = e (75*(x)-«7(x)) x ( x ) (2.21) 
f(x) = ^( x ) e (75 * {x)+ir?(x)) (2.22) 
V^W = J(<S>,r])VxVx (2.23) 

where $ and rj are scalar fields and J($, r/) is the Fujikawa Jacobian As it is well 
known, the above transformation permits to decouple the field A^ from the fermionic 
fields if one writes 

A M (x) = e^$(x) + <V/(x) (2.24) 

which can also be considered as a bosonic change of variables with trivial (field indepen- 
dent) Jacobian. Then, we find 

det(0 + 4 + 7oA) = J($, v) det(0 + 7 A) (2.25) 

After a suitable regularization the fermionic Jacobian in Fourier space (x — > p = 
(p ,Pi)) reads PU 

f 1 f d 2 p 

J{$,r}) = expj-y [-p?$(p)$(-p)+P?^(p)»7(-p)-2poPi$(p)^(-p 

(2.26) 
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The vacuum to vacuum functional is then expressed as 



Z = N J V^Vr] exp- 5 ^ 1 J($, rj) det{i$ + 7 A) 
where S[<p, rj] is obtained when one inserts (|2.24|) in S[A^] (see eq. (|2.20jl ): 



(2.27) 




p)-A(p) + 77(P)^(-P)^(P) + *(pM-p)C(p)] (2-28) 



with 



•A(P) 



B(P) 



C(P) 



^[po^) 1 (p)+Pi(^ { o) 1 (p) + ^) 

^[^(o) 1 (P)+^( y (i) 1 (P)-^ 
(^) 1 (P)-^) 1 (P)+^)- 



(2.29) 



We have then obtained a completely bosonized action for the non-equilibrium system 
which can now be used as the starting point to explore specific time-dependent interac- 
tions. 

3 Green's functions 

We are now ready to compute the Green's functions from eq. ([2.27)1 . by taking functional 
derivatives with respect to suitable sources. It is worthwhile noting that these Green's 
functions will be time ordered along the contour C, in the complex time plane. From 
now on we will restrict our study to contact electron-electron interactions in which the 
couplings Vo and V\ are constants in momentum space. 

Let us now focus our attention on the one-particle fermionic propagator: 



here the subscripts for the Green functions denotes the branches of the contour in the 
complex time plane on which the time coordinates reside. The subscripts + and — refer 




(3.1) 



where 




(3.2) 
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to fields defined in the upper and lower branches, respectively, corresponding to forward 
(+) and backward (— ) time evolution. 

Using the decoupling technique described in the previous section we can factorize the 
Green's function as a product of a purely fermionic Green's function and a bosonic part 
Be, exactly as it happens in the equilibrium case To be specific we shall consider the 
propagator corresponding to right-moving fermions (the left-moving component can be 
obtained in a completely similar way): 



Gc(x, t;x,t') = Gf c (x, t; x',i) Bg(x -x',t-t 



Rl 



(3.3) 



with 



B^(x—x ,t—t ) = exp <!. / d 2 p sin 



7T 



Po(t - t') +pi{x - x) 



g(p) - A(p) + iC(p) \ 

£(p) J ' 

(3-4) 

where A, B and C is given by eq.flOU and £ = C 2 - 4AB. Please note that we are 
now using a more explicit notation in the arguments of functions (x = (x,t)). Since we 
are dealing with perturbations that depend explicitly on space and time, up to now in 
an arbitrary way, the Green's functions are not expected to be translationally invariant 
neither in time nor in space. 

It is convenient to express B^(x,t) in the form: 



B*(x) = exp(/ 1 (x)+J 2 (x)) 



(3.5) 



where I\p are given by 

lb- a 



7T 2 2ab 
1 b - a + 1 

7T 

7T 2 ab 



J d 2 p sin 2 [-p. x] exp 



•A|j>i| 



ij>l+V 2 i)ij>l + v 2 pl) 



.1 



ipoPi 



<i 2 psin 2 [^p.x]e A|pi1 -^- ,, , 

2 (p 2 + p 2 )(p 2 + v 2 p 2 ) 



(3.6) 



(3.7) 



where we have introduce an ultraviolet cutoff A and two parameters a and b defined by: 

_1 = V(i) = (92 - 9a) 



a 
b- 1 



V, 



(0) 



72 + 94)- 



The integrals I\ 2 can be explicitly evaluated and one obtains 



A(x) 
/ 2 (x) 



2 



x 2 + (\t\ +Af 



A 2 



1 



-{K + K- X )\n 



x 2 + + A.f 



A 2 



-sgn(t) ( In 



A+ I t I — ix 
A + Itl + ix 



In 



A + f |t| — ix 
A + t>|t| + ix 



(3.8) 



where 



K = 



VT 7T y 1 + W 71 " + S-2/7T 



1 + ^)(1 - ^) = A /(l + ^ + ^)(l ^ ■ 



7T 7T V TT 7T 7T 7T 



(3.9) 



are the parameters usually defined in the field theory description of Luttinger liquids. 
In particular v is the velocity of the charge- density modes that dominate the low energy 
physics of these systems. 

Inserting eq. ()3.8|) in eq. (|3.5)l we finally obtain 

^ ^u\( A + t - ix ) A27 , a ^{A-t + ix) A 27 



B C (X, t) Q(t) ^ A + vt _ ^ ^ 2 + ^ + A ^ 2 ^ 7 +©( *)-^ A _ vt + ^ ^ 2 + ^ + 

(3.10) 

where 

1 =l(K + K- 1 -2) (3.11) 

Now we consider the fermionic piece Gj t c{x, t; x', t') that satisfies the Euclidean equa- 
tion: 

{id t + d x + A(x, t)}G ItC {x, t; x', f) = -5(x - x')5 c (t - f). (3.12) 

with 5 C given by (l2~T4l . 

For later convenience, from now on we shall work in Minkowski space (t — > it and 
G — > — G). At this point we also choose a particular form for the perturbation: 

A(x,t) = A sm{Qt)e(t)[e{x + a/2)-e(x-a/2)) (3.13) 

where G(t) is the Heaviside function. We then have a separable harmonic perturbation 
adiabatically switched on, consisting of a square barrier of width a which height A oscil- 
lates in time with frequency Q. A similar situation was recently considered by Komnik 
and Gogolin in their study of transport through a tunneling junction between a metallic 
lead and a Luttinger liquid [IT] . However two important differences should be stressed. 
First of all they did not include the effect of the adiabatic switching. Secondly, they did 
not specify a particular form for the spatial piece of the interaction, but they were able to 
expressed time-averaged quantities in terms of a parameter related to the strength of the 
interaction. As we shall see, by specifying the spatial dependence of the barrier we will 
be able to give explicit results for the TED, before a temporal mean value is considered. 
The Green's function Gi t c(x,t; x' ,£') satisfies: 

{d t + d x + \ sin(ftt)e(t) [&(x + a/2) - &(x - a/2)]}G I c (x, t; x' , t') = -5{x - x')8 c (t - 1'). 

(3.14) 
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Gi t c{x, t', x', t') can be factorized in terms of the free (translationally invariant) Green's 
function and a function to be determined that contains the effect of the broken transla- 
tional symmetry: 

G IiC (x, t; x', t') = G 0iC (x -x',t- t') exp[f3(x, t) - (3(x', f)\, (3.15) 
where Go t c(x — x r ,t — t') satisfies: 

{d t + d x }G 0iC (x -y,t-t) = -5(x - y)6 c (t - t). (3.16) 

The solution reads 

G oC (x-y t-t) = Q^*-*') \ | Q c(t'-t) 1 

2n [(x — y) — (t — t') + iao] 2n [(x — y) — (t — t') — iao] ' 

(3-17) 

where «o is a n infrared cut-off, and we have defined the Heaviside function on the contour 
in the usual way: 

e(t-t') t,t' on C + 

BAt-f)=\ n e(( '- t > y » C -, - (3.18) 
» t on C+, v on C_ v ' 

1 if on C+, t on C_ 

The function (3(x,t) obeys the following equation: 

{d +dj}P(x,t) = -A sin(m)e(t) [Q(x + a/2) - 6(ar - a/2)], (3.19) 
which solution is given by 

\ poo rco 

/5( x ,t) = — / dt'sm(nt) / dx'[Q(x +a/2)-Q(x' -a/2)}G (x-x',t-t), (3.20) 
2tt Jo J-oc 

where G is the inverse of the (d t + d x ) operator. 
We thus obtain 

P(x, t) = -i^Q(t){Q(-(x+a/2)) cos(tt(t-(x+a/2)))-Q(-(x-a/2)) cos(Q(t-(x-a/2)))+ 

+ cos(Qt)[Q(x + a/2) - Q{x - a/2)}}. (3.21) 

At this point we already have all the ingredients needed to build the Green's function 
of equation ()3.3|) . Indeed, using (|3.17|) and (j3.21|) in (|3.15j) we have the first factor in 
()3.3|) . which together with ()3.10|) give us the complete Green's function. But before we 
write down the final result two important observations are in order. First we note that 
both ()3.10j) and 1)3.17)1 are translationally invariant. The effect of the time dependent 
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interaction is completely contained in the exponential factors. It is then natural to 
reexpress Gq in the form 

G c (x,t;x',t') = Gl(x-x',t-t)exp[/3(x,t) - {3(x',t% (3.22) 

where G]y(x — x',t — t') is the product of (j3.10|) and (j3.17|) . which is nothing but the 
equilibrium propagator for a Luttinger liquid with its characteristic interaction dependent 
exponent 7. The second observation concerns the actual form of G^, which up to this 
point is the exact electronic propagator. As it is well-known, even in the equilibrium case, 
computations of physical quantities with this propagator are a non trivial challenge [Tlj . 
[To] . Therefore we shall consider the usual long-distance, and long-time approximation 
\x\, \t\ S> A which leads to 



Gl(x-x',t-t') = 

e c (t-t') r 1 



2n \ v{t — t') — (x — x') — ia 
<d c (t' - t) f 1 



A- 



(v(t - t') - (x - x') - iA)(v(t - t') + (x- x') - iA) 



2n \ v(t — t') — (x — y) + ia 



A 



2 



(v(t - t') - (x- x') + iA){v(t - V) + (x- x') + iA) 



(3 



Equation ()3.22j) . with the explicit forms for /5(x) and Gq given in (j3.21|) and ()3.23j) . 
is our first non trivial result. We have obtained an exact (within the long-distance 
approximation) analytical expression for the Green's function of a Luttinger liquid in the 
presence of a dynamic barrier. As we shall see in the next section, this result will enable 
us to compute the energy distribution of this system. 



4 Total Energy Distribution Function 

In this section we use the Wigner representation in which the distribution function is 
defined in terms of G_ + (x, x'). In this representation it is usually convenient to define 
a center of mass coordinate system (R,T) = (x + x')/2, (r, t) = (x — x'). One then 
Fourier transforms the relative variable, and in this way, in terms of the distribution 
function one obtains macroscopic quantities such as the particle density, particle current, 
energy density and energy current. Here we are specially interested in computing the 
total energy distribution (TED) of the electrons in the presence of the time-dependent 
perturbation defined in eq. (j3.13|) . The TED is defined as 

n(uj,x,t) = —i J dr exp(zc<jr) G_ + (x, t + r/2; x, t — r/2), (4.1) 

where we now call po = to. Of course, we can use the results obtained in the previous 
section for the Green's functions, in order to compute n(u, x, t) (see eq. ()3.2j) ). As 
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before, we shall consider the contribution to the TED of the right-moving particles only 
(ur(uj, x, t). When one tries to obtain an explicit expression for this quantity, it becomes 
apparent that the necessary stages of the calculation are greatly simplified if one considers 
the large-distance (time) approximation \x\, \t\ ^> A. As a first step, for illustrative 
purposes we will evaluate, within this approximation, the TED to first order in X/Q 
(computation of higher orders is tedious but straightforward). Thus we write 

n R (uj, x, t) = n (uj) - i^ni{u, x, t), (4.2) 
where the first term in the right hand side is the equilibrium TED given by 

nJu) = / dr- \\ - ' — (4.3) 

uv ; 2tt J {a + ivT)(A + ivT) 2 ~f v ; 

which, after taking the limit cto — > 0, yields 

. . exp(-Acu/v) ( Kuj\ 21 ^ ( AuA n , . , . 

where $(a, b; z) is Rummer's confluent hypergeometric function. This result agrees with 
previous calculations [Hj, In particular, in the low-frequency limit we reobtain the 
expected behaviour n (u) ~ uj 21 . 
The first non trivial correction is 

A 27 f°° dr 
n/(^,x,t) = — { / — — . — —[exp(i(uj + n/2)T)exp(int)F(x,a,tt) 

+ exp{i{uj — Q/2)r) exp(— iflt)F(x, a, — fl)] 
n2t dr 

-[exp(i(uj — VL/2)t) exp(iflt)F(x, a, Q) 



5 (a + ivr)(A + wt) 2 "< 
+ exp(i(uj + n/2)r)exp(-ittt)F(x,a,-tt)]} (4.5) 

which can be cast in the form 

A 27 f°° drexp(i(uj + n/2)T) 

4vr L v 1 v ' J_ 2t (a + ivr)(A + ivt) 2 ^ 

. \ t-i / ^ f°° drexpUiu -VL/2)t) 

+ exp(-iQt)F(x, a, -fl) \ r _ | , ,_ WA , + 

— exp(iQt)F(x, a, Q) 



2t (a + ivr)(A + zt>r) 2 T 
2 * drexp(i(cj - fi/2)r) 



oo 

2t 



a + wt)(A + wt) 2 ~< 



i dr expiiiuj + Vt / '2)t) , . A . 

- exp (-iQt )F (x, a, -fi) / '—^-f , 4.6 

' (a + zwr)(A + wr) 2 ')' J v 7 
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where 

F(x, a, fi) = 0(-(x + a/2)) exp[-itt(x + a/2)] - 6(-(x - a/2)) exp[-zfi(a; - a/2)] + 

+6(x + a/2) -0(x-a/2). (4.7) 

At this point we consider two important limits: for very large times in the past 
(t — > — oo) one obtains nj = 0, i.e. the TED recovers its equilibrium value, as expected. 
In the other case, for large values of t in the future, the integrals in the above expression 
coincide with ()4.3|) and one has: 

tii(uj, x, t) = -[exp(ifl,t)F(x, a, Q) — exp(— iQt)F(x, a, — fl)] [no(u + Q/2) — n^uo — Sl/2)\. 

2 (4-8) 
Then, up to first order in A, for t — > — oo 

n R (u,x,t) = n (u), (4.9) 

and for very large positive times: 

n R (co, x, t) = n {uo) — ^[exp(ifit)F(x, a, Q)— exp(— ittt)F(x, a, — Q)] [no(^+0/2)— n (uj— Sl/2)\. 

(4.10) 

Thus, in this case the non-equilibrium correction to the TED is a superposition of 
unperturbed TED's centered at energies ±Q/2 . 

In the large time approximation we were able to extend the above result to all orders 
in A. Indeed, taking into account that, in this regime one can write 

oo / '\\ n i 

exp(/?(x,t + r/2))-exp(/?(x,t-r/2)) « ^ f -L- J - (exp(^r/2) - exp(-^r/2)) n 

x (exp(zf2t)F(x, a, fi) — exp(— iQt)F(x, a, — f2)) n , 

and using the binomial expansion 

(exp(zfir/2) - exp(-zfir/2)) n = ^ ( ^ ) (-l) n ~ fc exp(iflr(" - 2fc)/2), (4.11) 

fc=0 ^ ' 

one finally gets 



00 n 1 / Ai \ n 

t) = V V — ( ^ (exp(ittt)F(x, a, SI) - exp(-ittt)F(x, a - tt)) ) x 
/ — ' n! \ 2iz / 

-l) n ^n (c; + fi(n - 2Jfe)/2)(4.12) 



n R (u;,x, 

n=0 fc=0 

n 
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Figure 1: : N(u, z = 0.1) for 7 = (filled line) and 7 = 1/2 (dashed line). 



which is one of the main results of this work. We derived, in the large time regime, an 
analytical expression for the TED at all orders in X/Q. 

In realistic systems the frequency Q is expected to be quite high so that it is unlikely 
that the explicit time resolution of the TED would be experimentally accessible. Then, 
it is natural to consider the time average, over the period of the perturbation, for our 
TED: < tir{oj,x) >= N(u,x). In so doing one readily discovers that only even powers 
of X/fl survive the mean value: 

= EE ( r, "J (L-k)M {{mj F ^a,n)F(x,a,-n)j n (u+n(n-k)) 

(4.13) 

Then, the averaged TED is a superposition of unperturbed TED's centered at energies 
nQ, with n integer. The terms in this linear combination are weighted by powers of 
z = (-^) 2 F(x,a,Q)F(x,a, —ft), i.e. factors associated to the strength of the interaction 
and its geometry. The interpretation of the physical consequences of the combined effect 
of the electron-electron interaction and the time dependent impurity is facilitated by 
evaluating N(u,x) in a particular spatial point and considering the function N(u, z), as 
was done in [TT]. It should be stressed, however, that our z-parameter differs from the 
one defined in that work, in which the absence of adiabatic switching allowed the authors 
to compute N(uj, z) without specifying a particular spatial dependence for the scatterer. 

In figures 1, 2 and 3, we show the behaviour of the function N(u, z) as function of 
u, for z = 0.1, 1 and 10, respectively. We have set O = v = A = 1, and recall that the 
Fermi energy has been put equal to zero. In each figure we have considered the effect 
of the dynamic barrier for both the free (7 = 0, filled line) and interacting (7 = 1/2, 
dashed line) cases. The results are in qualitative agreement with the ones obtained in jTT] 
for a generic time-dependent impurity without adiabatic switching. First of all one sees 
the appearance of equidistant sidebands in the TED corresponding to the noninteracting 
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Figure 2: : iV(u;, ,2 = 1) for 7 = (filled line) and 7 = 1/2 (dashed line). 
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Figure 3: : N(u, z = 10) for 7 = (filled line) and 7 = 1/2 (dashed line). 
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Figure 4: : N(x) for a = ir, 7 = (filled line) and 7 = 1/2 (dashed line). 
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Figure 5: : N(x) for a = 3n/2, 7 = (filled line) and 7 = 1/2 (dashed line). 



system |lfij . As it is known, this is a manifestation of the fact that, in the presence of a 
time dependent perturbation, scattering ceases to be elastic and gain or loss of energy can 
take place. When forward-scattering is taken into account the step functions of the free 
case become smoothed, and the values of N(u, z) are always below the ones corresponding 
to the noninteractine; expected. For increasing z, the number of sidebands also 

increases and the free and interacting values of N(u, z) become closer from each other. 

Since we have obtained not only the TED as function of u, but also its spatial de- 
pendence, we can also analyze the spatial distribution of energy for given energies. In 
figures 4, 5 and 6 we plotted N(x) for A = u — 1 and a = tt, 3tt/2 and 2rr respectively. 
As before, the filled line corresponds to the free case and the dashed line represents the 
behavior for 7 = 1/2. Due to the specific form of the impurity, the modulating factors in 
()4.13|) . which are determined by the function F(x,a, fi) (see equation ()4.7|1 ). vanish for 
x > a/2 (please remember that we are considering the time averaged TED for the right 
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Figure 6: : N(x) for a = 2ir, 7 = (filled line) and 7=1/2 (dashed line). 



movers only). For this reason N(x) acquires the unperturbed (equilibrium) value to the 
right of the barrier. On the other hand, the behavior of N(x) to the left of the scatterer 
is governed by the value of sin 2 (ft a/2). This is why when we vary a from n/Q to 2ii/fl 
the modulation decreases until it reaches the zero value, for a = 2n/VL. This last case 
corresponds to the symmetric distribution represented in figure 6, which is in fact the 
case of interest in realistic systems. Indeed, in carbon nanotubes [TZj and quantum wires 
[TH] , the width of the impurities are of the order of nms. Then, recovering the dependence 
on the Fermi velocity (which we had set equal to 1, but is of order 8 x 10 5 m/s) and 
considering external frequencies in the range of the far megahertzs, one concludes that 
Vta/{2v F ) w 10" 8 . 

Finally, in figures 7 and 8, we present three-dimensional plots of iV as function of 
both uj and x. 
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Figure 7: : N(lu, x) for a = 2n and 7 = 0. 




5 Conclusions 



We have undertaken the extension of a functional bosonization method, previously devel- 
oped in the context of equilibrium field theory, to the case in which time-dependent per- 
turbations are included. In particular, we have considered a model of fermions interacting 
through current-current interactions. In order to make contact with the description of 
ID condensed matter systems those interactions are incorporated in a non-covariant way. 
In addition to this, by now standard (equilibrium) setting, we also included a term in the 
action -bilinear in fermion fields- which describes scattering through a time-dependent 
barrier. By using the CTP (Schwinger-Keldysh) approach to manage the time evolu- 
tion of a field theory, we obtained the exact fermionic Green's function for our out of 
equilibrium Thirring model (equation (|3.22|) ). In the long-distance and long-time regime 
(the usual situation considered in the condensed matter context) we obtained a closed 
expression for the TED (equation (j4.12jl ). which is a superposition of equilibrium TED's 
centered at energies equal to integer multiples of the external half- frequency 0/2 . From 
this expression we readily obtained the time averaged TED N(u, x) (equation (j4.13|) ). in 
which only contributions corresponding to integer multiples of Q survive, modulated by 
factors related to the spatial distribution of the impurity. Our results are consistent with 
the findings of the authors of ref. who computed the temporal mean value of the 

TED for an arbitrary impurity geometry in a specific spatial point. In contrast to these 
authors, we included a temporal Heaviside function in the time dependent interaction. 
Although for this case we could not compute the TED for arbitrary scatterer geometries, 
by considering a rectangular barrier we were able to obtain explicit analytical results for 
the time-averaged TED, not only as function of energy but also in terms of the spatial 
coordinate. 

This work could be continued in several directions. Firstly, it would be desirable to 
include the effect of finite temperature. Secondly, one could try to explore the short- 
time regime by building the TED using directly equation ()3.22|) for the Green's function 
instead of the long-time approximation resulting through the simplified expression for 
Gq (equation ()3.23|) ). This does not seem to be an easy task, but nevertheless we hope 
to follow this route in the close future. 
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